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ABSTRACT 


-^■The  distribution,  abundance,  and  species  composition  o'" 
subtidal  macrobenthic  invertebrates  of  the  lower  Chesapeake 
Bay  ^-were  studied.  ^The  macrobenthic  infauna  was  sampled  four 
times  per  year  (seasonally)  at  12  sites  in  1982  and  1983. 
Sites  were  located  from  the  mouth  of  the  Chesapeake  Bay  along 
the  access  channels  to  the  port  of  Hampton  Roads,  within  the 
Elizabeth  River  and  up  the  James  River  to  Mulberry  Island. 
Samples  for  commercially  important  benthos  were  collected  at 
eight  locations  during  the  winter  and  summer  cruises  of  each 
year.  The  purpose  of  this  study  was:  (1)  to  present  recent 
information  concerning  the  structure  of  marcobenthic 
communities  of  the  lower  Chespeake  Bay,  (2)  to  compare  the 
data  generated  by  this  study  with  data  from  previous  studies 
of  the  lower  Chesapeake  Bay  and  estuaries  along  the 
Southeastern  U.S.,  and  (3)  to  develop  multivariate 
statistical  models  based  upon  the  baseline  data  and  test  the 
sensitivity  of  these  models  to  simulated  impacted  data  sets. 

j>From  a  cluster  analysis  of  the  infaunal  collection 
sites,  five  major  site  groups  were  defined.  The  12  infaunal 
collection  sites  were  organized  by  the  cluster  analysis  into 
groups  that  were  generally  spatially  contiguous.  This 
clustering  reflected  gradients  of  decreasing  average  grain 
size,  water  depth  and  salinity  that  occurred  moving  up  the 
estuary.  -  > 


''  Total  oomn  unity  density  and  indices  of  species 

■  r  ' 

di'/e»'sity  were  within  the  ranges  reported  in  previous  studies 
for  comparaPle  sedinent  types  and  salinity  ranges.  The  above 
oarameters  were  lowest  at  one  site  near  the  mouth  of  the  bay 
and  at  the  site  group  defined  by  two  sites  within  the 
E 1 izabeth  River. 

Qualitative  samples  for  commercially  important  benthic 
invertebrates  (blue  crabs,  hard  shell  clams  and  oysters) 
confirmed  d  i  s t r i b u t i ona 1  patterns  well  known  to  local 
f  i  s  h  e  rm  en.^a  nd  indicated  by  numerous  unpublished  surveys  of 
the  lower  bay. 

•Multivariate  models  were  developed  for  each  of  the  site 
groups.  The  sensitivity  of  each  type  of  model  was  tested 


using  simulated  impacted  data  sets 


INTRODUCTION 


The  distribution  and  abundance  of  the  macrobent  hie 
invertebrates  of  the  lower  Chesapeake  Say  were  studied. 
Density  dominants,  community  abundance,  species  diversity, 
and  animal-sediment  relationships  were  determined  from  data 
from  12  sites  extending  from  the  mouth  of  the  Chesapeake  Bay, 
along  the  access  channels  to  the  port  of  Hampton  Roads,  into 
the  Elizabeth  River  and  up  the  James  River  to  Mulberry 
Island.  These  sites  were  part  of  an  environmental  study 
concerning  the  potential  effects  of  deepening  the  access 
channels  of  the  lower  Chesapeake  Bay.  Samples  were  collected 
four  times  per  year  from  January  1982  through  October  1983 
for  the  macrobenthic  infaunal  community  and  twice  per  year 
for  commercially  important  benthic  species. 

The  purpose  of  this  study  was:  (1)  to  present  recent 
information  concerning  the  structure  of  macrobenthic 
communities  of  the  lower  Chesapeake  Bay,  (2)  to  compare  the 
data  generated  by  this  study  with  data  from  previous  studies 
of  the  lower  Chesapeake  Bay  (Boesch  1  97  3  ;  Dauer  et  al.  1984; 
Ewing  and  Dauer  1982;  Hawthorne  and  Dauer  1983;  Tourtellotte 
and  Dauer  1983)  and  other  estuaries  of  the  mid-Atlantic  and 
Southeastern  U.S.  (Dorjes  and  Howard  1975;  Maurer  1977; 
Maurer  et  al.  1978,  1979;  Tenore  1972),  and  (3)  to  develop 
multivariate  statistical  models  based  upon  the  baseline  data 
and  test  the  sensitivity  of  these  models  to  simulated 


impacted  data  sets. 


Materials  and  Methods 


The  macrofauna  of  the  lower  Chesapeake  Bay  was  sampled 
at  12  sites  located  from  the  mouth  of  the  bay,  along  the 
access  channels  to  the  port  of  Hampton  Roads,  through  the 
Southern  Branch  of  the  Eli2abeth  River  and  up  the  James  River 
to  Mulberry  Island  (Fig.  1).  Samples  were  collected  four 
times  per  year  (seasonally)  in  January,  April,  July  and 
October  of  1982  and  1983* 

All  samples  for  the  macrobenthic  infaunal  community  were 

2 

collected  using  a  Shipek  grab  (surface  area  of  0.04  m  ). 
During  the  initial  cruise  in  January  1903,  15  Shipek  grabs 

were  collected  from  sites  A,  C,  E,  G,  I  and  K  (Fig.  1)  in 
order  to  determine  the  sample  size  required  for  an  a  £r^or_i 
determined  level  of  precision.  Each  grab  was  washed  through 
a  0.5  mm  mesh-sized  screen,  relaxed  with  dilute  isopropyl 
alcohol,  and  preserved  and  stained  with  a  f o r m a  1 i n  -  r o se 
bengal  solution. 

The  number  of  Shipek  grabs  necessary  to  acquire  a 
statistically  reliable  estimate  of  the  density  of  individuals 
was  determined  using  the  following  formula: 

\Dx/ 

where:  s  =  standard  deviation  of  the  preliminary  sample, 

t  =  the  tabulated  t  value  at  the  0.05  level  with 
the  degrees  of  freedom  of  the  preliminary  set 
of  samples 


X  =  mean  density  of  the  preliminary  sample 

D  =  required  level  of  precision  expressed  as  a 
decimal  (Southward  1966) 

The  product  Dx  represents  the  half  width  of  a  desired 
95%  confidence  interval,  that  is,  the  distance  between  the 
actual  mean  and  the  field  measured  mean  with  a  probability  of 

0.95.  This  is  a  measure  of  the  precision  desired  in  our 

sampling  program.  Previous  work  with  benthic  organisms  has 
shown  that  a  distance  of  30  to  35  percent  of  the  mean  will 
give  a  statistically  reliable  estimate  (Dauer  et  al.  1979). 
With  a  30  percent  level  of  precision,  an  average  of  6.2 

Shipeks  per  site  would  be  necessary  (range  2.9  to  9.8)  and 
with  a  35  percent  level  of  precision  and  average  of  •’.8 

Shipeks  per  site  would  be  necessary  (range  2.2  to  7.5). 
Based  upon  these  calculation  and  considering  the  manpower 
available,  five  Shipeks  per  site  were  used  to  characterize 
the  benthic  infaunal  community. 

At  each  site  a  small  portion  of  the  surface  sediment 
(8  drams)  was  retained  for  sediment  analysis.  If  the 
sediment  from  an  individual  grab  changed  markedly,  an 
additional  sediment  sample  was  taken.  Dry  sieving  of  the 
sand  fraction  and  a  pipette  analysis  of  the  silt-clay 
fraction  were  conducted  using  the  techniques  of  Folk  (1974). 
Mean  particle  size,  sorting  coefficient  and  percent  sand  were 
determined  graphically  using  the  formulae  of  Folk  (1974). 

At  each  collection  site  on  each  cruise  bottom  salinity 
and  water  depth  were  determined. 


During  the  January  and  July  cruises  of  each  year 
commercially  important  benthic  species  were  sampled  from  10 
minute  tows  using  a  commercial  clam  dredge  at  3  sites  shown 
on  F i g  .  1  . 

Community  Ana -  dominant  Spec i^es 

All  infaunal  taxa  collected  were  used  in  the  computation 
of  indices  of  community  structure.  Shannon's  informational 
diversity  index,  Margalef's  species  richness  index,  and 
Pielou's  eveness  index  were  calculated  (see  Ewing  and  Dauer 
1982,  for  further  details). 

The  12  collection  sites  were  stratified  into 
biologically  determined  groups  based  upon  a  cluster  analysis. 
In  this  analysis  all  species  that  had  a  two  year  mean  density 
greater  than  0.5  individuals  per  grab  sample  over  all  sites 
were  included  in  the  analysis  (35  species).  Ta xono m i ca  1 1  y 
problematic  taxa  that  could  not  be  accurately  identified  to 
the  species  level  were  excluded  (e.g.  Oligochaeta  spp.  and 
Cirratulidae  spp).  The  selected  species  were  used  in  a 
normal  classification  analysis  of  the  sites  using  the  Bray- 
Curtis  similarity  coefficient  and  group  average  sorting  on 
logarithmically  transfomed  data  (Boesch  1977b). 

Because  some  species  were  collected  in  very  high 
numbers  in  only  a  few  sites  and/or  collection  times, 
dominance  of  species  for  the  entire  study  was  based  upon  the 
biological  index  ran'xing  (McCloskey  1970).  For  each  cruise 
the  top  ten  density  dominants  were  scored.  The  species  with 


the  highest  density  received  a  score  of  10,  the  species  with 
the  second  highest  density  received  a  score  of  9,  etc.  Rank 


density  scores  were  summed  over  all  cruises  and  all  sites 
within  a  site  group.  Only  the  top  10  species  were  used  in  the 
following  analyses. 

For  each  of  the  top  ten  density  dominants  of  each  site 
group  defined  by  the  cluster  analysis  the  minimum  amount  of 
change  necessary  to  produce  a  statistically  significant 
difference  was  determined  using  the  following  formula  (from 
Sokal  and  Rohlf  1969  as  modified  by  Michael  et  al.  198l); 


[^2s^(t^(v)  ^2  (1  _p  )  (  V  )  )  ^  ^ 


where: 


G  =  smallest  detectable  true  difference 
3  =  standard  deviation 

(  V  )  ’ 2  (  1 -P  )  (  V  )  '  values  from  a  two-tailed  t  table  with  v 

degress  of  freedom  and  probabilities  of 
■x-  and  2(1-P),  respectively 
P  =  desired  probability  that  the  difference 
will  be  detected 

tX  =  desired  level  of  statisitlcal  significance 


Multivariate  Models 

Statistical  models  were  developed  and  the  "sensitivity" 


m  i 

*?V 


•  i 


O  O 


V.  '  •  *  -  >.  ^  •  *-  M*.  -%  _S  n"  Vi  "-iK  '  •  \ 


7 


(i.e.,  the  ability  of  the  model  to  detect  changes  of  known 
values)  of  these  models  to  a  variety  of  potential  impacts  was 
tested  using  simulated  impacted  data  sets  (SIDS).  For 
further  discussion  of  the  rationale  of  this  approach  see 
Alden  et  al.  (1982). 

The  SIDS  were  produced  by  a  computer  program  developed 
by  Dr.  R.W.  Alden  III.  For  each  species  used  in  the  analysis 
the  SIDS  were  generated  to  have  the  same  frequency 
distribution  as  the  baseline  data,  but  with  different  mean 
values  that  represented  potential  impacts.  Briefly,  the 
program  used  a  power  law  transformation  to  produce  the  best 
fit  to  the  baseline  data,  changed  the  true  mean  to  a  desired 
mean,  and  then  untransformed  the  data.  Any  desired  number  of 
replicates  could  be  produced.  SIDS  were  produced  with  the 
same  auto-  and  crosscorrelation  relationships  as  the  baseline 
data.  For  each  of  the  site  groups  defined  by  the  cluster 
analysis  6  different  SIDS  were  produced  as  follows:  each 
species  reduced  in  mean  density  by  50{,  75%,  or  90%;  each 
species  increased  in  mean  density  by  50%,  100%,  or  200%.  For 
each  type  of  impact  five  impacted  data  sets  were  produced. 

Stepwise  discriminant  analysis  was  used  to  develop 
models  to  test  for  differences  between  groups  defined  a 
pr^orj^.  Two  groups  were  defined  -  one  group  was  one  of  the 
cluster  analysis  site  groups  while  the  second  group  was  one 
of  the  SIDS.  Discriminant  analysis  produces  a  multivariate 
linear  additive  model  that  best  d e s c r i m i na t e s  between  the 


defined  groups.  The  model  is  then  tested  by  classifying  all 


replicates  (baseline  and  SIDS)  into  one  of  the  two  groups, 
and  checking  the  percentage  of  correct  classifications.  The 


optimal  model  will  classify  all  replicates  from  the  baseline 
data  into  one  group  and  all  replicates  of  the  SIDS  into  the 
other  group  (lOOj  correct  classification).  For  sensitivity 
testing  a  significant  impact  was  declared  to  have  occurred  if 
greater  than  95%  of  the  SIDS  replicates  were  correctly 
classified . 

A  second  type  of  model  was  based  upon  the  approach 
suggested  by  Green  (1979)  for  baseline  monitoring  studies.  A 
principal  components  analysis  was  conducted  upon  each  of  the 
site  groups  defined  by  the  cluster  analysis.  A  principal 
components  analysis  produces  a  multivariate  linear  additive 
model  with  the  first  principal  component  accounting  for  the 
greatest  amount  of  variance  in  the  data  set.  The  next 
principal  component  is  independent  of  the  previous  one  and 
accounts  for  the  greatest  amount  of  residual  variance.  This 
process  is  continued  for  all  remaining  principal  components. 
Green's  approach  produces  a  two-dimensional  graph  based  upon 
the  first  two  principal  components.  A  95%  probability  elipse 
Is  calculated  for  standardized  principal  component  scores  for 
the  first  two  principal  components.  SIDS  are  next  compared 
to  the  baseline  data  and  the  difference  in  principal 
components  scores  are  plotted.  If  a  plotted  point  lies 
outside  the  probabilty  ellipse  a  significant  impact  is 
indicated;  if  within  the  probability  ellipse  no  impact  is 
ind ica  ted  . 
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Results 


S^te  Characteristics 

The  study  sites  ranged  from  moderately  sorted  medium 
sands  (Site  D)  to  very  poorly  sorted  fine  silts  (Site 
J)(Table  1).  As  expected  sediments  tended  to  become  finer 
in  grain  size  and  more  poorly  sorted  moving  up  into  the 
estuary.  Salinity  and  water  depth  also  tended  to  decrease 
moving  up  the  estuary  (Table  1). 

Community  Analyst s  -  Dominant  Species 

A  total  of  227  taxa  were  identified.  Polychaetes 
comprised  45. (101  species)  of  the  fauna,  bivalves  12, St 
(28  species),  amphipods  11. 5t  (26  species)  and  gastropods 
11. Ot  (25  species).  See  Appendix  for  a  complete  listing  of 
the  taxa  identified. 

From  the  cluster  analysis  five  site  groups  were  defined 
(Fig.  2).  Table  2  shows  the  major  sedimentary  and  physical- 
chemical  data  as  averaged  over  all  sites  within  a  site  group. 
In  general,  proceeding  from  site  group  1  to  5  sediments 
become  finer  in  grain  size  and  more  poorly  sorted  while 
salinity  and  water  depth  decrease.  This  pattern  reflects 
primarily  the  spatial  arrangement  of  the  site  groups  with 
site  group  1  near  the  mouth  of  the  bay  and  site  groups  4  and 
5  farthest  up  the  estuary.  Site  groups  1-3  form  a  "sand" 
grouping  with  mean  particle  values  in  the  fine  and  very  fine 
sand  categories.  Site  groups  4  and  5  form  a  "mud"  grouping 
with  mean  particle  values  in  the  coarse  to  fine  silt 


categories  (Table  2). 

Table  3A  summarizes  total  community  parameters  as 
calculated  by  the  major  site  groups.  Diversity  and  total 
community  density  were  lowest  in  site  group  5  which  consisted 
of  the  two  sites  in  the  Elizabeth  River. 

The  density  dominants  for  the  five  site  groups  are 
shown  in  Tables  4-8. 

^mmer c^aj^  Ben tho£ 

Table  9  lists  the  commercially  important  benthic 
invertebrates  collected  from  the  eight  dredge  sites  (Fig.  1). 
The  results  confirm  well  known  patterns  of  distribution  for 
commercial  species  of  the  lower  bay  as  indicated  by  numerous 
unpublished  studies  and  information  from  local  fisherman. 

The  blue  crab  ( Ca 1 1 Inec te s  ia£idu£)  was  collected  in 
the  open  bay  collection  sites  (2  and  3)  almost  exclusively  in 
winter  samples  (96. 0?  of  individuals  in  winter  samples  only). 
The  open  bay  collections  represent  samples  of  overwintering 
females  that  form  the  basis  for  the  crab  dredging  industry  of 
the  lower  bay.  The  other  large  collections  of  blue  crabs 
were  in  the  Elizabeth  River  sites  (6  and  7,  Fig.  1)  with 
individuals  collected  primarily  (93-2%)  in  summer  samples. 

The  hard  shell  clam  ( Mercenar ia  mercenarla)  was  highly 
concentrated  in  dredge  sample  sites  3-5  (Table  9)  again 
reflecting  known  distributional  patterns  for  commercially 
Important  populations  of  the  species  (Mayne  et  al.  1982). 

The  highly  concentrated  distribution  of  the  oyster 
®  dredge  site  8  corresponds  with 


collections  of  "seed  bed"  populations  of  the  James  River 
(Table  9). 

Several  species  of  minor  commercial  importance  were 
also  collected  (Table  9).  The  rock  crab  (Cancer  irroratus) 
and  blue  mussel  ( My tilus  e  d  u 1 i s )  are  much  more  commercially 
important  farther  north  along  our  coast  and  do  not  currently 
represent  important  resources  for  the  Chesapeake  Bay. 


Multivariate  Models  and  Sensitivity  Analysis 

Table  10  shows  the  results  of  the  sensitivity  analysis 
using  the  discriminant  models.  For  a  given  simulated  impact 
the  mean  percent  correct  classification  is  shown.  For  all 
site  groups  an  increase  of  100-2001  in  all  species  would  be 
necessary  to  declare  a  significant  impact,  while  only  site 
groups  3  and  4  were  very  sensitive  to  large  decreases  in 
density. 

For  all  site  groups  a  total  defaunation  was  never 
declared  to  be  significant  (i.e.  was  always  plotted  within 
the  95%  probability  ellipse).  Therefore,  the  principal 
components  model  proved  to  be  insensitive,  and  unacceptable 
for  future  impact  assessment.  Increases  in  the  range  of 
500-800<  were  necessary  to  indicate  a  significant  impact. 

Tables  3-8  also  show  the  minimal  detectable  impact  (as 
a  percent  of  the  mean  value)  for  total  community  parameters 
and  the  dominant  species  for  each  site  group.  For  the  total 
community  parameters  the  three  species  diversity  indices  are 
the  most  sensitive  parameters  with  total  community  density 


being  the  least  sensitive  (Table  3)»  Site  groups  1  and  2 
would  require  over  lOOt  change  in  the  mean  value  before  a 
significant  difference  could  be  detected.  This  basically 
reflects  the  higher  seasonal  variability  at  these  site  groups 
rather  than  between  replicate  variation  within  a  season. 

An  examination  of  patterns  of  M.D.I.'s  between  the  site 
groups  shows  that  site  groups  1  and  5  are  the  least  sensitive 
(or  most  variable)  statistically  to  changes  in  mean  values  of 
the  dominant  species  (Tables  4-8).  The  density  dominants  of 
these  two  site  groups  must  change  on  the  average  168. 3X  and 
102. 8t,  respectively,  be  fore  a  statistically  significant 
difference  could  be  detected.  Site  group  1  contains  a  single 
species  and  site  group  5  four  species  with  M.D.I.'s  less  than 
loot.  Site  group  4  is  the  most  sensitive  with  eight  species 
with  M.D.I.'s  less  than  100?  (Table  7)  and  an  average  M.D.I. 
of  80.61.  Site  groups  2  and  3  are  somewhat  intermediate  in 
sensitivity  with  average  M.D.I.'s  of  dominant  species  of 
93.91  and  92.31,  respectively. 


DISCUSSION 


Comparison  wi^  other  studies 

Dauer  et  al.  (1984)  compared  recent  studies  of  the 
raacrobenthos  of  the  lower  bay  (Dauer  et  al.  1984,  Ewing  and 
Dauer  1982;  Hawthorne  and  Dauer  1983;  Tourtellotte  and  Dauer 
1983)  to  studies  based  on  data  collected  approximately  10 
years  prior  (Boesch  1973*  1977a).  Comparisons  with  prior 
studies  were  made  difficult  by  the  fact  that  most  of  the 
previous  studies  in  the  Chesapeake  Bay  used  a  1.0  mm  screen 
in  comparison  with  the  0.5  mm  screen  used  in  the  more  recent 
studies.  Most  of  the  species  that  were  dominants  only  in  the 
more  recent  studies  ( S  trebl_03£io  bened^c_tl.,  Med  ioma  s  t  us 

gaperatus,  Te^^^na  ajj^^i^)  -^re  mostly 
smaller  sized  species  (most  or  all  individuals  pass  through  a 
1.0  mm  screen)  whose  abundance  estimates  would  be  greatly 
affected  by  the  screen  size  used  in  collections. 

For  the  five  site  groups  defined  in  this  study,  total 
community  density  and  indices  of  species  diversity,  species 
richness  and  eveness  were  generally  high  and  within  the 
ranges  reported  for  comparable  sediment  types  and  salinity 
ranges  (Table  3  of  this  study  compared  to  studies  reviewed  in 
Dauer  et  al.  1  984).  The  lower  densities  and  indices  of 
species  diversity  found  within  the  Elizabeth  River  site  group 
of  this  study  are  within  the  ranges  reported  in  previous 
studies  (Hawthorne  and  Dauer  1983)  and  generally  reflect  the 
lower  diversity  and  density  expected  in  high  silt-clay  type 


sediments  associated  with  creeks  and  rivers  of  limited 
circulation  (Dauer  et  al.  198^). 

There  are  several  species  previously  reported  as 
dominants  that  are  rare  or  absent  in  the  more  recent  studies 
-  Am£e^^sca  vadorum,  Heteromas^us  f f ££!!ii3 ,  and  Ne£ht^s 
sandy  habitats;  S££0 c ha e£o££e r u s 
£££21 'S££*lii£  and  L^no£herus  amb^gua  in  muddy 
sediments  (see  Dauer  et  al.  1984  for  further  comparisons). 
Reduction  in  the  dominance  of  these  species  may  be  due  to 
long  term  trends  of  the  various  species  or  to  major  episodic 
events.  Salinity  reduction  due  to  Tropical  Storm  Agnes 
produced  great  differences  in  density  of  the  macrobenthio 
dominants  of  the  polyhaline  region  of  the  James  River  and 
Hampton  Roads  area  (Boesch  et  al.  1976). 

The  dominant  species  of  this  study  (Tables  4-8)  can  be 
divided  into  four  groups:  Group  I  -  species  dominant  on  the 
inner  shelf  and  sandy  sediments  of  the  lower  bay,  Group  II  - 
species  dominant  in  sandy  sediments  of  the  lower  bay,  Group 
III  -  species  dominant  in  sandy  and  muddy  sediments,  and 
Group  IV  -  species  dominant  in  muddy  sites  only. 

Group  I  species  were  dominants  both  on  inner  shelf  and 
sandy  sediments  of  the  lower  bay  and  included  S£^o£hane£ 
bombyx,  Te 1 1 ina  agilis,  Nephtys  picta  and  Polygordlus  sp.  S. 
bombyx  and  agilis  are  the  most  widespread  species  of  this 
group  being  reported  as  dominants  of  the  inner  shelf  (Table 
7,  Dauer  et  al.  1984),  a  variety  of  locations  in  sandy 
sediments  of  the  open  bay  (bay-wide  transects,  clean-sand 


site  group,  Table  3,  Dauer  et  al.  1984;  eastern  shore, 
offshore  site  group,  Table  3,  Ewing  and  Dauer,  1982; 
Lynnhaven  Roads  site  group,  Table  6,  Tourtellotte  and  Dauer, 
1983)  and  associated  with  inlet-shoal  sites  of  lower  bay 
tidal  creeks  and  bays  (Table  4,  Ewing  and  Dauer,  1982;  Table 
3,  Tourtellotte  and  Dauer  1983).  ^  was  reported  as  a 
dominant  in  all  three  sand  groups  of  this  study,  while 
Polygordius  sp.  was  a  dominant  in  site  groups  1  and  2  of  this 
study • 

Group  II  species  were  dominant  in  sandy  sediments  of 
the  lower  bay  only  and  included  EnsJ.^  Aniasti^os 
££££.£££.££  and  Theses  species  were 
dominants  in  site  groups  1  and  2  (but  not  3)  of  this  study 
and  were  also  reported  as  dominants  of  the  clean-sand  site 
group  of  the  open  bay  study  of  Dauer  et  al.  (1984).  G. 
d ibranchiata  was  also  reported  as  a  dominant  species  of  sand 
sites  off  the  eastern  shore  of  the  bay  (Table  3,  Ewing  and 
Dauer  1982). 

Group  III  species  were  found  in  both  sandy  and  muddy 
site  groups  and  included  Mediomastus  ambiseta,  Paraprlonospio 
plnnata,  Glycinde  solitarla  and  Acteocina  canallculata.  M. 
ambiseta  is  the  most  widely  distributed  species  in  the  lower 
bay.  This  species  has  been  reported  as  a  density  dominant  on 
the  shelf  (Dauer  et  al.  1984),  of  all  site  groups  except  1  of 
this  study,  in  numerous  locations  and  sediment  types  of  the 
open  bay  (Tables  2  and  4,  clean-sand  and  silty-sand  site 
groups,  Dauer  et  al.  1984;  Table  3,  offshore  site  group, 
Ewing  and  Dauer  1982;  Lynnhaven  Roads  site  group.  Table  6, 


Tourtellotte  and  Dauer  1933)  and  in  muddy  sites  of  tidal 
creeks  and  small  coastal  bays  (Table  5,  Ewing  and  Dauer  1982; 
Table  5,  Tourtellotte  and  Dauer  1983).  Both  ^  canaliculata 
and  Pj_  P i.n n a ^  were  previously  described  by  Boesoh  (1973)  as 
ubiquitous  with  respect  to  sediment  type  in  the  Hampton  Roads 
area.  Both  species  were  reported  as  dominants  in  open  bay 
silty-sand  sites  (Table  4,  Dauer  et  al.  1984)  and  P.  pinnata 
as  a  dominant  in  high  silt-clay  content  sites  of  the 
Lynnhaven  (Table  5,  Tourtellotte  and  Dauer  1983)  and 
Elizabeth  Rivers  (Table  4,  Hawthorne  and  Dauer  1983).  ^ 

was  a  dominant  species  in  the  Lynnhaven  River 
system  in  offshore  fine  sands  and  in  muddy  creek  bottoms 
(Tables  5  and  6,  Tourtellotte  and  Dauer  1983). 

Group  IV  species  were  dominant  only  in  the  high  silt- 
clay  sediments  of  site  groups  4  and  5  of  this  study  and 
include  Leucon  am  er^  can  us,  S^rebj^0  3£io  benedict^,  Nereis 
££,££££  heteropoda  and  Leitoscoloplos  AH 

of  the  polychaetes  in  this  group  have  been  termed  "euryhaline 
opportunists"  by  Boesch  (1977a).  All  five  of  the  species  in 
this  group  were  previously  reported  as  dominants  of  the 
Southern  Branch  of  the  Elizabeth  River  (Table  4,  Hawthorne 
and  Dauer  1983).  ^  benedict!,  H.  succ inea  and  ^  heteropoda 

were  previously  reported  as  dominants  of  muddy  sediments  of 
the  eastern  shore  and  Lynnhaven  River  system  (Table  5,  Ewing 
and  Dauer  1982;  Table  5,  Tourtellotte  and  Dauer  1983). 

Comparisons  with  other  estuaries  of  the  mid-Atlantic 
and  Southeastern  U.S.A.  show  broad  qualitative  similarities. 
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Boesch's  (1977a)  estuarine  opportunisits  are  widely 
distributed  along  the  coast,  particularly  in  silty  sediments 
(Watling  197U;  Dcirjes  and  Howard  1975;  Maurer  1977;  Tenore 
1972).  The  fauna  of  sandy  habitats  of  the  polyhaline  regions 
of  Delaware  bay  (Maurer  et  al.  1978,  1979),  the  Pamlico  River 
(Tenore  1972)  and  Georgia's  Ogeechee  River  (Dorjes  and  Howard 
1  975  )  often  are  characterized  by  such  species  as  S£i^o£hane3 
bombyx,  Glycera  dibranchiata  and  Tellina  spp.  Ens^s  directus 
is  also  a  sand  species  of  wide  geographic  distribution 
(Maurer  et  al.  1  9  7  8  ,  1  9  79  ).  ^^Iieis  succ^nea  is  a 
characteristic  species  of  mesohaline  regions  (Diaz  and  Boesch 
1977;  Dorjes  and  Howard  1975;  Tenore  1972).  Macoma  balthica 
was  also  reported  as  a  dominant  species  of  the  mesohaline 
region  of  the  James  River  of  the  Chesapeake  Bay  (Diaz  and 
Boesch  1977)  and  in  the  Pamlico  River  (Tenore  1972). 

3nd  Sensitivity  Anal^^s^£ 

The  testing  of  the  sensitivity  of  various  multivariate 
models  developed  from  the  baseline  data  is  useful  (1)  to 
indicate  the  magnitude  of  change  necessary  to  produce  a 
statistically  significant  difference  and  (2)  to  test  if 
models  may  be  relatively  insensitive,  and  therefore, 
inappropriate  for  impact  assessment. 

The  discriminant  models  produced  for  the  different  site 
groups  varied  in  sensitivity  to  a  variety  of  simulated 
impacts  (Table  10).  For  this  study  95%  or  better  correct 
classification  of  the  simulated  impacted  data  sets  was 
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considered  to  indicate  a  significant  impact.  All  five  site 
groups  were  significantly  affected  by  a  2001  increase  in  the 
density  of  dominant  species,  while  only  site  group  3  was 
sensitive  to  a  100?  increase.  With  decreasing  densities  only 
site  groups  3  and  4  were  sensitive  and  only  when  densities 
were  decreased  at  least  901t. 

The  sensitivity  patterns  shown  are  probably  affected  by 
the  size  of  the  site  group,  and  therefore,  the  number  of 
replicates  used  to  create  the  discriminant  models.  Site 
groups  1  and  5  which  were  generally  the  least  sensitive 
consisted  of  one  and  two  sites,  respectively,  while  site 
groups  3,  4  and  5  consisted  of  three  sites  each  (Fig.  2'. 

Although  site  groups  3  and  4  showed  the  greatest 
overall  sensitivity  to  simulated  impacts,  this  should  not  be 
confused  with  the  biological  sensitivity  of  the  component 
species  and  communities.  A  greater  absolute  degradation  of 
the  environment  at  site  groups  3  and  4  may  be  necessary  to 
produce  the  same  relative  amount  of  change  in  the  biota 
compared  to  site  groups  that  are  less  sensitive  to 
statistical  changes. 

The  graphical  method  of  Green  (1979),  which  is  based 
upon  a  principal  components  model  of  the  baseline  data,  was 
shown  to  be  too  insensitive  to  be  useful  in  impact 
assessment.  A  principal  components  analysis  produces  models 
useful  for  indicating  which  potential  factor(3)  might  explain 
or  account  for  the  greatest  amount  of  variance  in  the  data 
set.  However,  the  two-dimensional  graphical  technique  of 
Green  does  not  declare  a  total  defaunation  as  being 


significant:  such  a  result  is  ecologically  unacceptable. 
However,  without  the  type  of  sensitivity  testing  used  in  this 
study  this  model  may  have  been  used  in  future  impact 
assessment  studies.  In  that  case  ecologically  unacceptable 
alterations  might  occur  which  the  model  would  say  were  not 
statistically  significant.  The  necessity  of  sensitivity 
testing  with  simulated  impacted  data  sets  is  obvious. 


SUMMARY 


1.  The  m a c ro b en t h i c  infaunal  communities  associated 
with  access  channels  to  the  port  of  Hampton  roads  and  the 
James  and  Elizabeth  Rivers  were  sampled  quantitatively  in 
January,  April,  July  and  October  of  1982  and  1983  at  12 
sites.  Commercially  important  benthic  invertebrate  species 
were  sampled  at  eight  sites  during  the  January  and  July 
cruises  of  each  year. 

2.  From  a  cluster  analysis  of  the  infaunal  collection 
sites,  five  major  site  groups  were  defined.  Each  site  group 
was  analyzed  separately.  The  original  collection  sites  were 
organized  by  the  cluster  analysis  into  groups  that  were 
generally  spatially  contiguous.  These  cluster  site  groups 
reflected  gradients  of  decreasing  average  grain  size,  water 
depth  and  salinity  progressing  from  the  mouth  of  the  bay  up 
the  estuary. 

3.  Total  community  density  and  indices  of  species 
diversity,  species  richness  and  eveness  were  generally  high 
and  within  the  ranges  reported  in  previous  studies  of  the 
lower  bay  for  comparable  sediment  types  and  salinity  ranges. 
The  lowest  values  for  these  parameters  were  found  at  one  site 
near  the  mouth  of  the  bay  and  at  the  site  group  composed  of 
two  sites  within  the  Elizabeth  River. 

U.  The  dominant  species  of  all  site  groups  were 
classified  into  four  species  groups  based  upon  their 
distributions.  Group  I  species  were  dominants  on  both  the 
inner  continental  shelf  and  sandy  sediments  near  the  mouth  of 


the  bay.  Group  II  species  were  dominants  in  sandy  sediments 
near  the  mouth  of  the  bay.  Group  III  species  were  dominant 
in  both  sandy  and  muddy  sediments  (ubiquitous  species). 
Group  IV  species  were  dominant  only  in  muddy  sediments. 

5.  Qualitative  samples  for  commercially  important 
benthic  invertebrate  species  confirmed  well  known 
distributional  patterns  for  blue  crabs,  hard  shell  clams  and 
oysters . 

6.  Multivariate  models  based  upon  discriminant  and 
principal  components  analysis  were  developed  for  each  site 
group.  The  sensitivity  of  each  model  was  tested  using 
simulated  impacted  data  sets. 

7.  The  minimal  amount  of  change  necessary  in  order  to 
detect  a  significant  difference  was  calculated  for  all 
community  parameters  (total  community  density  and  species 
diversity  indices)  and  for  each  density  dominant  of  each  site 


group . 
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Table  1 


Summary  of  sedimentary  and  phys'cal-chemical  data  for 
the  study  sites.  A.  Sec  mentary  particle 
distribution  parameters.  B.  Salinity  (in  ppt)  and 
water  depth.  Shown  are  means  for  all  parameters. 


A.  Sedimentary  Parameters 


SITE 


MEAN  PHI 

SORTING 

COEFFICIENT 

t  SAND 

2.44 

0.87 

95.47 

3.  17 

1.29 

65 . 48 

2.63 

0.90 

90.04 

1 . 80 

0.80 

94.33 

3.33 

0.74 

86.54 

4.21 

1.50 

37.60 

3.08 

0.95 

77.82 

4.67 

1  .49 

21.70 

6.62 

3.05 

12.27 

7.17 

2.98 

6.19 

5.67 

2.79 

14.17 

4.53 

1 .08 

14.78 

B.  Salinity  and  Water  Depth 


SITE 


SALINITY  (ppt) 

29 . 1 

29.7 

28. 1 
28.0 

25.1 
24.5 

20.7 

19.1 


DEPTH  (m) 


Table  2. 

Summary  of 
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and 

physical 

-chemical  data 

by 
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site  groups  defined  by 

the 

cluster 

analysis. 
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Table  9.  Summary  of  the  spatial  distribution  of  commercially 
important  benthic  species.  Shown  are  the  total 
numbers  of  individuals  collected  in  four  10  minute 
trawls  of  a  clam  dredge  (January  and  July  cruises  of 
1982  and  1983).  A.  Major  commercial  species. 

B.  Minor  commercial  species. 


A.  Major  Commercial  Species 

Collection  Site 


Species 

1 

2 

3 

4 

5  6 

7 

Callinectes  sapidus 

- 

13 

12 

1 

5  13 

31 

Mercenaria  mercenaria 

- 

- 

25 

6 

33 

- 

Crassostrea  virginica 

- 

- 

- 

- 

3 

- 

B.  Minor  Commercial  Species 

Collec  tion 

Site 

Species 

1 

2 

3 

4 

5  6 

7 

Cancer  irroratus 

8 

7 

4 

- 

- 

- 

My tl lus  edulis 


100  750  300 


Table  TO.  Summary  of  sensitivity  analysis  of  discriminant  models 


based  upon  baseline  data  and  simulated  impacted  data  sets. 
Sliown  are  the  average  percent  correct  classification  of  the 
simulated  impacted  data  sets.  Each  percent  shown  is  the 
average  of  5  simulated  impacted  data  sets. 


Simulated  Impact 

1 

2 

3 

4 

5 

2005 

i nc  rea  se 

96.0 

97.3 

98.7 

97.3 

100.0 

1005 

increase 

80.0 

75.3 

96.0 

90.7 

74 . 0 

505 

i nc  rea  se 

80.0 

64 . 0 

74 . 7 

80 . 0 

72.0 

505 

decrease 

84 . 0 

85.3 

93.3 

81.3 

86 . 0 

755 

decrease 

80 . 0 

86.7 

92.0 

85.3 

88 . 0 

905 

decrease 

88.8 

84.0 

94.7 

94.7 

88.0 

APPENDIX 


TAXA  COLLECTED  DURING  STUDY 


CNIDARIA  ;  ANTHOZOA 

Anthozoa  spp. 

PLATYHELMINTHES  :  TURBELLARIA 
Turbellaria  spp. 

NEMERTEA 

Nemertea  spp. 

ANNELIDA  :  POLYCHAETA 

Aglaophamus  clrclnata  (Verrill) 
Amastlgos  oaperatus  Ewing  and  Dauer 
Ampharete  arotioa  Malmgren 
Arapharetldae  spp. 

Anclstrosyl 11s  hartmanae  Pettibone 
Apoprionosplo  pygmaea  (Hartman) 
Arlcldea  Gather Inae  Laubier 
Art Gldea  wassi  Pettibone 
Asabellides  oculata  (Webster) 

Asyohis  e longata  (Verrill) 

Autoly tus  spp. 

Branla  clavata  (Claparede) 

Brania  welf leetensls  Pettibone 
Capitella  oapltata  (Fabricius) 
Capltella  spp . 

Capltomastus  aolculatus  Hartman 
Cirratulidae  spp. 

Clrrophorus  furoatus  (Hartman) 
Clymenella  torguata  (Leidy) 

Dlapatra  ouprea  (Bose) 

Drllonerels  longa  Webster 
Drllonerels  magna  Webster  and  Benedict 
Eteone  heteropoda  Hartman 
Eteone  lactea  Claparede 
Eumlda  sangulnea  (Oersted) 

Glycera  amerlcana  Leidy 
Glycera  dibranohiata  Ehlers 
Glycera  spp. 

Glycinde  sol i tarla  (Webster) 

Gyptls  brevipalpa  ( Har tmann-Schroder ) 
Gyptls  vlttata  Webster  and  Benedict 
Harmothoe  extenuata  (Grube) 

Hemlpodus  roseus  Quatrefages 
Heteromastus  f iliformls  (Claparede) 
Hydroldes  d lanthus  Verrill 
Leltoseoloplos  f ragills  (Verrill) 
Lepldonotus  aublevls  Verrill 
Lolmla  medusa  (Savlgny) 

Lumbrinerls  f ragi lis  (Muller) 

Lumbrlner is  tenuis  Verrill 
Macroclymene  zonalls  (Verrill) 

Magelona  sp. 

Maldanidae  spp. 

Medlomastus  amblseta  (Hartman) 


t 
t 

Micropthalmus  slmilis  Sobretsky 
Mlnusplo  cirri  fen  ( W  i  r  e  n  ) 

Nephtyidae  spp. 

Nephtys  inclsa  Malmgren 
Nephtys  picta  Eh 1  era 
Nereidae  spp. 

Nereis  acuminata  E h 1 e r s 
Nereis  suocinea  (Frey  and  Leuckart) 
Notomaatus  hemipodua  Hartman 
Notomastus  laterlceus  Sars 
Onuphidae  spp. 

Onuphia  ereml ta  Audouin  and  Milne-Edwarda 
Owenla  fuaiformis  delli  Chiaje 
Paleanotua  heteroseta  Hartman 
Paradoneia  ly ra  (Southern) 

Paraheaione  luteola  (Webster) 

Paranaitia  apeciosa  (Webster) 

Parapr lonoaplo  £in n a t a  ( E h 1 e r a  ) 

Pectlnaria  gouldii  (Verrill) 

Pheruaa  a p . 

Phyllodoce  arenae  Webster 
Phvllodoce  caatanea  (  Marenze Her ) 


ANNELIDA 


ANNELIDA 


MOLLUSCA 


MOLLUSCA 


Streblospio  b e n e d Webster 
Sy llldae  spp . 

Syllides  ve rr 1 1 1 1  Moore 
Terebellidae  spp. 

:  OLIGOCHAETA 
Oligochaeta  spp. 

;  HIRUDINEA 
Hirudinea  spp. 

:  gastropoda 

Acteocina  canallculata  (Say) 

Anachi s  obesa  Adams 
Busycon  oarica  (Mont  fort) 
Corambella  depressa  B a  1 c h 
Cory phe 11a  s p . 

Crepldula  fornlcata  ( L i n n e  ) 
Cyclostremiscua  ^auiA  (  P i s o h e r  ) 
Cyllchnella  bldentata  ( 0 r b 1 g n y ) 
Epltonlum  multlstrlatum  (Say) 
Epltonlum  ruplcola  (Kurtz) 
Epltonlum  s p . 

Eup leura  caudata  (Say) 

Gastropoda  spp. 

Mange  11a  cer Ina  Kurtz  and  Stlmpson 
Mltrella  lunata  (Say) 

Nassarlus  trl vlttatus  (Say) 
Nassarlus  vlbex  Say 
Natlca  pusllla  Say 
Odostomla  spp. 

Pollnloes  duplicatus  (Say) 

Polycera  s p . 

Rlctaxls  punctostrlatus  (Adams) 
Turbonl 1 la  Interrupta  (Totten) 
Turbonl 1 la  spp. 

Turrldae  spp. 

;  BIVALVIA 

Allgena  elevata  (Stlmpson) 

Anadara  ovalls  (Brugulere) 

Anadara  transversa  (Say) 

Blvalvla  spp. 

Cerastoderma  plnnulatum  (Conrad) 
Chlone  cancellata  Linnaeus 
Crassostrea  vlrginlca  (Graelln) 
Ensls  d Irectus  Conrad 
Eucrassatel la  speclosa  (Adams) 
Gemma  gemma  (Totten) 
lachad lum  recurvum  (Raflnesque) 
Lyonsla  hyallna  Conrad 
Macoma  balthlca  Linnaeus 
Macoma  mltchelll  Dali 
Macoma  tenta  Say 
Mercenarla  mercenarla  (Llnne) 

Mul Inla  lateralis  (Say) 

Mya  arenarla  Linnaeus 

My sel la  planulata  (Stlmpson) 

Mytllus  edulls  Llnne 


Nucula  proxima  Say 
Pandora  bushiana  Dali 
Pandora  trlllneata  Say 

Parvilucinia  multillneata  (Tuomey  and  Holmes) 

Siliqua  costata  Say 

Spisula  solldlsslma  ( D i 1 1 w y n ) 

Telllna  agl 1 1 s  Stirapson 
Yoldia  limatula  (Say) 

ARTHROPODA  :  ISOPODA 

Cyathura  pol 1 ta  (St imp son) 

Edotea  triloba  (Say) 

Erlchsonella  flllformls  (Say) 

Ptl lanthura  tenuis  (Harger) 

ARTHROPODA  :  AMPHIPODA~ 

Ampellsca  abdita  Mills 
Ampelisca  vadorum  Mills 
Ampellsca  verrilll  Mills 
Batea  catharinensls  Muller 
Caprellldae  spp. 

Cerapus  tubularis  Say 
Corophlum  spp. 

Elasmopus  levi s  Smith 
Erlchthonius  brasiliensls  (Dane) 

Gammarus  mucrona tus  Say 
Leptooheirus  plumulosus  Shoemaker 
Llstriella  barnardl  Wigley 
Llstriella  olymenellae  Mills 
Mellta  appendiculata  (Say) 

Melita  nl tlda  Smith 
Parametopella  cyprls  (Holmes) 

Parapleustes  aestuar lus  Watling  and  Maurer 
Photls  macrocoxa  Shoemaker 
Pleustidae  sp. 

Protohaustorlus  spp. 

Stenothoe  mlnuta  Holmes 
Stenothoe  s p . 

Synchelldlum  amer Icanum  Bous field 
Trlchophoxus  eplstomus  (Shoemaker) 

Unclola  Irrorata  Say 
Unclola  serrata  Shoemaker 
ARTHROPODA  :  CUMACEA 

Cyclaspls  varlans  Caiman 
Leucon  amerlcanus  Zimmer 
Oxyurostylls  smi thi  Caiman 
Pseud oleptocuma  minor  (Caiman) 

ARTHROPODA  :  MYSIDACEA 

Mysidopsis  blgelowl  Tattersall 
Neomysis  amer Icana  (Smith) 

ARTHROPODA  :  DECAPODA 

Callinectes  sapidus  Rathbun 
Cancer  irroratus  Say 
Crangon  septemspinosa  Say 
Decapoda  spp. 

Euceramus  praelongus  Stimpson 

Hexapanopeus  angustif rons  Benedict  and  Rathbun 


Liblnia  emarginata  Leach 
Ogyrides  limlcola  Williams 
Ovalipes  ooellatus  ( H e r b s t ) 
Pagurus  s p p . 

Panopeus  herbstil  Milne-Edwards 
Plnnixa  chaetopterana  Stimpson 
Pinnixa  cristata  Rath bun 
Plnnixa  sayana  Stimpson 
Pinnotheridae  spp. 

Upogebla  af finis  (Say) 

Xanthidae  spp. 

ARTHROPODA  :  PYCNOGONIDA 
Pycnogonida  spp. 

SIPUNCULA 

Phascolion  s tromb i  (Montagu) 

ECHIURA 

Echiura  spp. 

PRIAPULIDA 

Priapulida  spp. 

PHORONIDA 

Phoronls  psammophl la  Cori 
ECHINODERMATA  ;  ASTEROIDEA 

Aster ias  forbesil  ( D e s o r  ) 
ECHINODERMATA  :  ECHINOIDEA 

Arbacia  punctulata  (Lamarck) 
Echinarachnius  parma  ( L a r m a c k ) 
ECHINODERMATA  :'"hOLOTHUROIDE  A 
Holothuroidea  spp. 

Leptosynapta  inhaerens  (Ayres) 
ECHINODERMATA  ;  OPHIUROTdEA 
Ophiuroidea  spp. 

HEMICHORDATA 


